HIGH-MASS CLOUD CORES IN THE r] CARINAE GIANT 

MOLECULAR CLOUD 



YOSHINORI Y0NEKURA\ ShIN 'ICHIRO ASAYAMA^'^, KiMIHIRO KlMURA^, 

Hideo Ogawa\ Yoko Kanai^, Nobuyuki Yamaguchi^'^, 
Yasuo Fukui^, and Peter J. Barnes^'^ 

O 
O 

(N 

ABSTRACT 

■ We carried out an unbiased survey for massive dense cores in the giant molecular cloud 

'\0 ! associated with r/ Carinae with the NANTEN telescope in ^^CO, "CO, and C^^O J = 1-0 

emission lines. We identified 15 C^*0 cores, whose typical line width AT4omp, radius r, mass 
' M, column density -/V(H2), and average number density n(H2) were 3.3 km s~^, 2.2 pc, 2.6x10'^ 

Mq, 1.3x10^^ cm~^, and 1.2x10^ cm~^, respectively. Two of the 15 cores are associated with 
I IRAS point sources whose luminosities are larger than 10"*^ Lq, which indicates that massive star 

formation is occuring within these cores. Five cores including the two with IRAS sources are 
OO ; associated with MSX point sources. We detected Hi3C0+ (J = 1-0) emission toward 4 C^^O 

cores, two of which are associated with IRAS and MSX point sources, another one is associated 
only with an MSX point source, and the other is associated with neither IRAS nor MSX point 
. sources. The core with neither IRAS nor MSX point sources shows the presence of a bipolar 

' molecular outflow in ^^CO (J = 2-1), which indicates that star formation is also occuring in the 

core, and the other three of the four H^^CO^ detections show wing-like emission. In total, six 
O I C^*0 cores out of 15 (= 40%) are experienced star formation, and at least 2 of 15 (= 13 %) 

^ are massive-star forming cores in the 77 Car CMC. We found that massive star formation occurs 
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^ I preferentially in cores with larger A^(Il2), M, 7i,(H2), and smaller ratio of M^-„/M. We also 

found that the cores in the 77 Car CMC are characterized by large AV and M^i^/M on average 
compared to the cores in other GMCs observed with the same telescope. These properties of the 
^ . cores may account for the fact that as much as 60-87 % of the cores do not show any signs of 

^ ' massive star formation. We investigated the origin of a large amount of turbulence in the 77 Car 

CMC. We found that turbulence injection from stellar winds, molecular outflows, and supernova 
remnants which originated from stars formed within the CMC, are not enough to explain the 
existing turbulence. We propose the possibility that the large turbulence was pre-existing when 
the CMC was formed, and is now dissipating. Mechanisms such as multiple supernova explosions 
in the Carina flare supershell may have contributed to form a CMC with a large amount of 
turbulence. 

Subject headings: ISM: clouds — ISM: individual(r; Carinae CMC) — ISM: molecules — radio lines: ISM 
— stars: formation 
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(GMCs), while low- mass stars form both in dark 
clouds and GMCs. The formation process of an 
isolated low-mass star is fairly well understood; 
the gravitational collapse of a molecular cloud into 
a protostellar core and the subsequent accretion 
via a circumstellar disk surrounded by an infalling 
envelope (e.g., Shu, Adams, & Lizano 1987 and 
references therein), followed by the dissipation of 
momentum into its surroundings through molecu- 
lar outflows (Lada 1985; Fukui et al. 1986, 1993; 
Fukui 1989; and references therein). The current 
paradigm is now being tested by observations, and 
it has become clear that some modifications to the 
details of protostellar formation are needed (e.g., 
Andre et al. 1993; Onishi et al. 1996, 1998; Evans 
1999 and references therein). 

Concerning massive star formation, an empir- 
ical scenario of the later stages of the massive 
star formation process has been presented (e.g., 
Kurtz 2000 and references therein); dense (~ 10^ 
cm~'^), compact (< 0.1 pc), and hot (> 100 K) re- 
gions called Hot Cores are formed in GMCs, then 
the UV radiation from newly formed stars within 
Hot Cores ionizes the surrounding material to pro- 
duce ultra-compact H ii regions (UCHHs), and 
finally the central stars become visible, such as 
the Trapezium cluster in Orion. On the contrary, 
evolution before the formation of Hot Cores re- 
mains unclear, although some attempts to detect 
sources at the earliest stages of massive star forma- 
tion have been made; surveys toward IRAS point 
sources and maser sources (e.g., Beuther et al. 
2002a; Mueller et al. 2002; Sridharan et al. 2002; 
Shirley et al. 2003) have revealed the presence of 
the dense, compact, and massive cores (density 
> 10^ cm^^, temperature > 100 K, luminosity > 
10** Lq) sometimes called High-Mass Protostellar 
Objects (HMPOs; Sridharan et al. 2002), some of 
which are not accompanied by UCHHs. Although 
sources without UCHHs may be younger than the 
rest, the presence of IRAS point sources and maser 
sources themselves stongly indicates that stars are 
already formed within HMPOs. It is quite natural 
to suppose that dense, compact, massive, starless 
cores (density ~ 10^ cm~^, size < 1.0 pc, mass > 
10'^ Mq) are the precursors of the HMPOs and Hot 
Cores (e.g., Evans et al. 2002). Systematic studies 
of these massive starless cores had not been carried 
out until recently, due to the difficulty in select- 
ing targets from the existing database. It is to 



be noted that large-scale surveys for possible can- 
didates for massive starless cores (density ~ 10'' 
cm~^) have already been made; e.g., Lada, Bally, 
& Stark (1991) surveyed the L1630 cloud (Orion 
B) in CS (J = 2-1). The importance of the study 
of starless cores has been revealed in a series of 
unbiased surveys for dense cores in low-mass star- 
forming regions, especially from a viewpoint of set- 
ting observational constraints on the evolutional 
time scales of the theoretical model (Onishi et al. 
1996, 1998; Tachihara et al. 2002). The same is 
expected for the massive star forming regions and 
thus unbiased surveys for massive dense cores in 
massive star-forming regions are also crucial. 

One of the most important tasks in the study of 
massive star formation at the moment is to iden- 
tify massive dense cores which have not yet ex- 
perienced massive-star formation but are capable 
of future massive-star formation. High-resolution 
observations of these cores will give observational 
constraints on the initial phases of the evolu- 
tional model of a massive star; massive stars are 
formed by direct accretion onto a central proto- 
star through an accretion disk (e.g., McKee & 
Tan 2003), which is simple extension of low-mass 
star formation, or by collisions of intermediate- 
mass protostars within a young cluster (e.g., Bon- 
nell. Bate, & Zinnecker 1998). Moreover, such 
observations are the only method to reveal the 
initial conditions of massive star forming cores, 
since the properties of the cores change rapidly 
soon after the formation of massive stars due to 
strong UV radiation and stellar winds. Surveys 
for Hot Cores and HMPOs toward IRAS point 
sources and maser sources are, however, biased to- 
ward regions which have experienced star forma- 
tion, and thus it is not possible to obtain informa- 
tion on the initial conditions of the massive star 
forming cores. Therefore studies of dense cores 
which have not yet experienced massive-star for- 
mation are awaited. Unbiased surveys for massive 
dense cores over entire GMCs in molecular lines 
as well as in the far-infrared and (sub-) millimeter 
continuum may be the most effective methods to 
search for massive starless cores which are capable 
of future massive-star formation. Far-IR and sub- 
millimeter/millimeter continuum observations are 
best suited to detect radiation from the dust in 
starless cores, since the radiation from cold dust 
(T < 20 K) peaks at far-IR (> 150 fxm) wave- 
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lengths with gradual decrease toward longer wave- 
lengths. Recent development of sensitive instru- 
ments, such as bolometer array receivers, makes 
(sub-) millimeter continuum observations more ef- 
fective for carrying out a large-scale survey for 
massive dense cores. An unbiased survey of en- 
tire GMCs with (sub-) millimeter continuum has 
not yet been made, but several surveys in small 
fields surrounding IRAS point sources are now 
in progress for the purpose of detecting massive 
dense cores located near the IRAS point sources. 
So far, regions of ~ 10' x 10' around ~ 200 
IRAS point sources have been observed at a reso- 
lution of 8-24" (Forbrich et al. 2004; Faiindez et 
al. 2004). Based on these surveys, a few massive 
dense cold cores without mid-IR (MSX) and far- 
IR (IRAS) counterparts have been found (Forbrich 
et al. 2004; Garay et al. 2004). Molecular line ob- 
servations have the invaluable merit of drawing 
kinematic information within the cores, which is 
impossible in continuum observations. The empir- 
ical relation between the luminosity of the forming 
star and the line width of the parent dense core has 
been found by many authors from observations in 
NII3 (e.g., Stacy et al. 1988; Wouterloot, Walms- 
ley, & Henkel 1988; Myers, Ladd, & Fuller 1991; 
Harju, Walmsley, & Wouterloot 1993; Ladd, My- 
ers, & Goodman 1994; Jijina, Myers, & Adams 
1999) as well as in C^^O (e.g., Saito et al. 2001), 
which makes molecular line observations still more 
important. 

We, the NANTEN group, are now conducting 
unbiased surveys for dense cores in the C^^O (J 
= 1-0) line with the NANTEN 4-meter millime- 
ter telescope. Tachihara et al. (2002) compiled 
a sample of 174 C^^O dense cores in low-mass 
star-forming regions of Taurus (Onishi et al. 1996, 
1998), L1333 (Obayashi et al. 1998), Ophiuchus 
North (Nozawa et al. 1991), p Oph (Tachihara, 
Mizuno, & Fukui 2000), Chamaclcon (Mizuno et 
al. 1999), Lupus (Hara et al. 1999), Corona Aus- 
tralis (Yonekura et al. 1999), the Pipe nebula (On- 
ishi et al. 1999), and the Southern Coalsack (Kato 
et al. 1999) including intermediate-mass cluster- 
forming regions such as p Oph, Cha I, R CrA, 
and Lupus 3, and found that cores with active 
star formation tend to have larger column den- 
sity, number density, and mass, and tend to be 
gravitationally more bound. Observations toward 
massive star forming regions are also in progress; 



we have detected 112 cores so far in the unbiased 
surveys for dense cores in GMCs such as Orion 
A, Orion B, Cepheus 0B3, Vela, S35/S37, and 
Centaurus (Yu, Nagahama, & Fukui 1996; Naga- 
hama 1997; Yamaguchi et al. 1999a; Saito et al. 
1999, 2001; Aoyama et al. 2001). The beam size 
of the NANTEN telescope, 2'.7 (-0.3 pc at the 
distance of the nearest GMC, Orion Ji, D = 400 
pc), is not necessarily sufficient to resolve massive 
starless cores. However, it is still meaningful to 
select candidates for massive starless cores (den- 
sity ~ 10^ cm~^, size ~ a few pc, mass > a few 
xlO'' Mq, which arc estimated from the physical 
parameters of massive starless cores), which can 
be traced by C^^O (J = 1-0) observations with 
NANTEN, since only a few massive starless cores 
have been identified so far. In this paper, we re- 
port the results of an unbiased survey for dense 
cores in the GMC associated with the rj Carinae 
nebula (the 77 Car GMC hereafter), with the aim 
of obtaining samples of massive- stoHes5-core can- 
didates which are possible sites for future massive 
star formation. 

The T] Car GMC has been one of the most active 
star-forming regions in the Galaxy at least until 
the formation of the youngest stellar cluster asso- 
ciated with the GMC (~3 Myrs ago, e.g., Feinstein 
1995). Eight stellar clusters (Trumpler 14, 15, 16, 
CoUinder 228, Bochum 10, 11, NGC 3293, and 
NGC 3324; hereafter we abbreviate Trumpler 14 
and Trumpler 16 as Tr 14 and Tr 16, respectively) 
are associated with the GMC, which contains more 
than 64 0-type stars including one of the most 
massive stars, ry Car (Feinstein 1995). When we 
estimate the star formation activity from the rich- 
ness of the clusters by counting the number of the 
most massive stars (e.g., stars with the spectral 
type of 03) since the total number of the clus- 
ter members can not be counted, the 77 Car GMC 
(with five 03 stars) is estimated to be the second 
most active star-forming region after NGC 3603 
with six 03 stars (Maiz-Apellaniz et al. 2004). Its 
proximity to the Smi {D — 2.5 kpc for the rj Car 
GMC compared to 6.9 kpc for NGC 3603) makes 
the T] Car GMC the best site for the detailed study 
of massive star formation, because much higher 
spatial resolution can be achieved. There is a pos- 
sibility that active star-forming regions, such as 
W49, where more than 100 0-type stars have been 
found in infrared observations (Alves & Home- 
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ier 2003), arc hidden in the galactic plane due to 
heavy foreground extinction. However, it is nat- 
ural to think that regions with heavy foreground 
extinction arc located distant from the Sun (11.4 
kpc in the case of W49), and thus the t] Car GMC 
still remains advantageous. It is to be noted that 
the well-studied massive star-forming region, the 
Orion region, is a much smaller system, contain- 
ing only a single 04-6 star (Maiz-Apellaniz et al. 
2004). 

Several observations have been made so far cov- 
ering the entire 7] Car GMC with moderate spatial 
resolution. Grabelsky et al. (1987, 1988) observed 
the whole ri Car GMC (3x2 deg^) in the ^^CO 
(J = 1-0) line with an 8'.8 beam, and revealed 
the existence of a string of CO clouds between I 
= 284?7 and 289° with a total mass of 6.7 x 10^ 
Mq. Zhang et al. (2001b) made sub-millimeter 
observations in ^^CO (J = 4-3) and atomic car- 
bon [CI] ^Pi-^Po with a 3' beam, and found that 
these two emissions and ^^CO (J = 1-0) are ap- 
proximately coextensive, and that the excitation 
temperature is higher (~ 50 K near rj Car and ~ 20 
K in clouds farther from r] Car) than that in dark 
clouds. High-resolution observations at a resolu- 
tion of ~20-40" were made toward some selected 
clouds near r] Car. There are two main CO emis- 
sion regions near rj Car; northern and southern 
clouds (de Graauw et al. 1981; Brooks, Whiteoak, 
& Storey 1998). The northern cloud lies surround- 
ing the cluster Tr 14 (~ 7 pc northwest of r] Car). 
Brooks et al. (2003) observed the northern cloud 
and found that the cloud consists of 4 main CO 
clouds with masses in the range 40-500 Mq, two 
of which are gravitationally bound with density 
exceeding 10^ cm~^. They proposed that these 
clouds arc exposed to a radiation field originating 
mainly from the most massive stars in Tr 14. The 
southern cloud is located adjacent to the cluster 
Tr 16, which includes rj Car. Megeath et al. (1996) 
found a gravitationally-bound cloud with a mass 
of 70 Mq, and that a cold IRAS point source with 
luminosity ~ 10"* Lq is associated with the cloud. 
They concluded that star formation is on-going at 
the interface between the cloud and the cluster Tr 
16, which is triggered by radiation driven shocks. 
In between the northern and the southern clouds, 
the distribution of the molecular gas is patchy due 
to dissociation by the strong UV radiation and 
stellar winds from nearby massive stars in Tr 16 



(Cox & Bronfman 1995). They also found that the 
typical mass of the clouds is ~ 10 Mq, and that the 
clouds are not gravitationally bound. Recently, 
Rathborne et al. (2004) studied the region called 
the giant piUar (Smith et al. 2000), ~ 0?5 south 
of T] Car. Since the distance from the clusters 
is larger than that of the northern and southern 
clouds, and the giant pillar is located behind the 
southren cloud, the influence of the UV radiation 
and stellar winds from the clusters may be smaller 
and the photo-dissociation takes longer than the 
northern and southern clouds. Prom CO J = 2- 
1 observations toward four mid-IR sources identi- 
fied by Smith et al. (2000) using mid-infrared data 
obtained by Midcourse Space Experiment (MSX), 
they found these sources are also influenced by the 
radiation from clusters Tr 14 and Tr 16. The mass 
of the associated clouds ranges from 10-45 Mq 
and they are probably not gravitationally bound. 
They also found 12 candidates for massive young 
stellar objects and UCHIIs from the mid-IR data. 
These studies revealed that star formation is still 
on-going in the 77 Car GMC under the strong in- 
fluence of the most massive stars in clusters, in at 
least 1.5 X 2 deg^ regions in the vicinity of 77 Car. 

Here we report the results of an unbiased sur- 
vey for dense cores in the entire 77 Car GMC. The 
purpose of the present study is to obtain a sample 
of massive- starless-core candidates which are pos- 
sible sites for future massive star formation. With 
moderate spatial resolution, 2.'7, our observation 
covered the entire r] Car GMC (3x2 deg^) in- 
cluding regions farther from r] Car, such as the 
western part of the rj Car GMC around the H 11 
region Gum 31 and the southern most part of the 
7] Car GMC, where no previous observations had 
been made except coarse large-scale surveys. Ob- 
servations were made in C^^O (J = 1-0) as well 
as in i^co and ^^CO with NANTEN. We also 
made high-resolution observations in H^'^CO^ (J 
= 1-0) with the Mopra 22-meter millimeter tele- 
scope and searched for molecular outflows in ^^CO 
(J = 2-1) toward some selected dense cores with 
the ASTE 10-meter sub-millimeter telescope. Out 
of the 15 C^^O cores detected in this survey, two 
arc on-going massive star forming cores associated 
with luminous infrared sources and another 4 are 
possibly less-massive star-forming cores, and thus 
as many as 9-13 cores which have not experienced 
massive star formation, are detected in this sur- 
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vey. We present the physical properties of the 
dense cores and compare the star-formation ac- 
tivity with them. We also compare the proper- 
ties of the C^*^0 cores with those in other mas- 
sive star forming regions, such as Orion, Cepheus, 
Vela, S35, and Centaurus, in order to study the 
characteristics of the t] Car GMC. 

2. OBSERVATIONS 

2.1. Large Scale Observations of CO Gcis 
and the Search for High-Mass Cores 

Observations of the J = 10 transitions of 
i^CO, i^CO, and C^^^O were made with the NAN- 
TEN 4-meter telescope of Nagoya University at 
Las Campanas Observatory of the Carnegie Insti- 
tution of Washington (Fukui et al. 1991; Fukui & 
Sakakibara 1992; Fukui & Yonekura 1998). The 
half-power beam width of the telescope was 2f7 at 
110 GHz and 2'.6 at 115 GHz, corresponding to 
~2 pc for the distance of the rj Car GMC, 2.5 kpc. 
The 4 K cooled SIS mixer receiver provided a typ- 
ical system temperature of ~140 K (SSB) at 110 
GHz and ~200 K (SSB) at 115 GHz, including 
the atmosphere toward the zenith (Ogawa et al. 
1990). The spectrometer was an acousto-optical 
spectrometer with a total bandwidth of 40 MHz 
divided into 2048 channels. The effective spectral 
resolution was 40 kHz, corresponding to a velocity 
resolution of 0.11 km s~^ at 110 GHz. 

The 12 CO data were obtained from 1998 May 
to July, and in 1999 March. Most of the region 
was observed by using a position-switching tech- 
nique with a grid spacing of 2', and the rest with 
4' grid. The total integration time and the typical 
rms noise of the data at a velocity resolution of 
0.1 km s~i were ~5 sec and AT^ms ^ 1-0 K, re- 
spectively. In -"^"^CO, observations were made from 
1999 January to March with 4' grid by using a 
position-switching technique. The total integra- 
tion time per point was ~5 sec and the typical 
rms noise of the data was ATrms ~ 0.4 K at a ve- 
locity resolution of 0.1 km s~^. C^^O observations 
were made in 3 periods from 1998 July to August, 
in 1999 March, and from 2003 August to October. 
The data were obtained with a grid spacing of 2'. 
A frequency-switching technique was used with a 
switching interval of 13 MHz in 1998 and 1999, 
whereas a position-switching technique was used 
in 2003. The total integration time per point was 



~90 sec and the typical rms noise of the data was 
ATrms ~ 0.1 K at a velocity resolution of 0.1 km 

For the calibration of the spectral line inten- 
sity, a room-temperature chopper wheel was em- 
ployed. The absolute-intensity calibration was 
made by observing Ori KL (aigso = 5'^32"47!0, 
5i950=-5°24'21") and p Oph East (aigso = 
16^29'"20!9, (5i95o = -24°22'13"). The peak ra- 
diation temperatures, T^, are assumed to be 
(i^CO, Ori KL) = 65 K (Ulich & Haas 1976; Kut- 
ner & Ulich 1981; Levreault 1988), ("CO, Ori 
KL) 10 K (Mizuno et al. 1995), and (C^^O, 
p Oph East) = 4.4 K. The last value is consis- 
tent with the assumption that {C^^O, TMC-1 
[ai95o = 4h38'"42!0, (^igso = 25°35'45"]) = 2.0 K 
(Onishi et al. 1996), which has been established 
by the long-term monitoring of standard sources 
made at the NANTEN telescope. The pointing 
accuracy was measured to be better than 20", as 
checked by optical observations of stars with a 
CCD camera attached to the telescope as well as 
by radio observations of Jupiter, Venus, and the 
edge of the Sun. 

2.2. Detailed Observations for High-Density 
Regions 

Observations of the J = 1-0 transition of 
R^^CO+ were made with MOPRA, the 22-meter 
radio telescope of the Australia Telescope Na- 
tional Facility (ATNF), on 2003 July. The half- 
power beam width of the telescope was ^37" at 
86 GHz (Ladd et al. 2005). The front end was a 
dual-channel 4 K cooled SIS mixer receiver, and 
the typical system temperature was 200 K (SSB) 
at 86 GHz, including the atmosphere toward the 
zenith. 

Two orthogonal polarizations were simultane- 
ously observed, both tuned to the E.^^CO~^ (J = 
1-0) frequency. The spectrometer was a digital 
correlator with a bandwidth of 64 MHz and 2048 
channels, which was divided into two 1024 chan- 
nels in order to measure the two polarizations si- 
multaneously. The velocity coverage was ~ 110 
km s~i at 86 GHz, and the effective spectral reso- 
lution was ~ 0.13 km at 86 GHz, correspond- 
ing to ~ 1.2 times the channel increment. 

The data were obtained with a grid spacing of 
40" by the position-switching mode. The total in- 



5 



tegration time per point was ^ 8 minutes, and 
the typical rms noise of the data was ^Trms 
0.1 K at a velocity resolution of 0.1 km s~^. The 
intensity calibration was made by using a room- 
temperature chopper-wheel. The absolute inten- 
sity calibration was made by observing Ori KL, by 
assuming = 0.55 K (Johansson et al. 1984). 

The pointing accuracy was measured to be 
within 15" by observing the SiO maser source 
RW Vel at 86 GHz every two hours during the 
observations. 

2.3. Search for Molecular Outflows 

Observations of the J = 2-1 transition of ^^CO 
were made with Atacama Submillimeter Telescope 
Experiment (ASTE), the 10-meter radio telescope 
of Nobeyama Radio Observatory'^ at Pampa la 
Bola, Chile, on 2003 November (Kohno et al. 2003; 
Ezawa et al. 2004). The half-power beamwidth of 
the telescope was 30" at 230 GHz. The front end 
was a 4 K cooled SIS mixer receiver (Sekimoto et 
al. 2001). The typical system temperature was 300 
K (DSB) at 230 GHz, including the atmosphere to- 
ward the zenith. We used a digital correlator with 
a bandwidth of 128 MHz and 1024 channels (Sorai 
et al. 2000). The effective spectral resolution was 
150 kHz, corresponding to a velocity resolution of 
0.16 km s"^ at 230 GHz. 

The data were obtained with a grid spacing of 
30" by using the position-switching mode. The to- 
tal integration time per point was ~ 1 minute and 
the typical rms noise of the data 
0.8 K at a velocity resolution of 0.5 km s~^. The 
intensity calibration was made by using a room- 
temperature chopper-wheel. The absolute inten- 
sity calibration was made by observing Ori KL by 
assuming T,^ to be 111 K (Sutton et al. 1985). The 
pointing accuracy was measured to be within 15" 
as checked by optical observations of stars with a 
CCD camera attached to the telescope as well as 
by radio observations of Jupiter and Mars. 

Nobeyama Radio Observatory (NRO) is a branch of the Na- 
tional Astronomical Observatory of Japan (NAO J) , which 
belongs to the National Institutes of Natural Sciences 
(NINS) 



3. RESULTS 

3.1. Distributions of CO Gas 

The J — 1-0 emission lines of carbon monoxide 
(i^CO) and its isotopes (i^CO and C^^O) have 
been used to study molecular clouds, since they 
are among the most abundant molecules in the 
dense interstellar medium. Although all those 
lines have nearly the same critical density for col- 
lisional excitation, it is found empirically that 
the typical densities traced by them arc different: 

- 102 cm--* for i^CO, - lO'"* cm--' for ^^CO, and 

- lO'' cm--"* for C^^^O (e.g., Mizuno et al. 1995, 
Onishi et al. 1996, 1998 for Taurus; Mizuno ct al. 
1998, 1999 for Chamaeleon; Dobashi et al. 1994, 
1996, Yonekura et al. 1997, Kawamura et al. 1998, 
Kim et al. 2004 for a series of ^^CO survey in the 
Milky Way; Tachihara et al. 2002 and references 
therein for C^^O cores in nearby star forming re- 
gions). This density differentiation can be under- 
stood as follows: One of the causes for this density 
differentiation is due to the optical depth. The 
^^CO transition is optically thick in most molec- 
ular clouds, making it to probe the lower den- 
sity enveloping layer of a molecular cloud due to 
the photon trapping. The other two lines have 
smaller optical depths, allowing us to look more 
deeply into clouds. However, because of the se- 
lective photo-dissosiation of the ^^CO and C^^O 
molecules, abundance ratio of the isotopes vary 
as a function of the visual extinction Ay (e.g., 
Warin, Benayoun, & Viala 1996). The abundance 
ratio R(^^CO/C^^O) reaches up to ~ 10 near the 
cloud surface at Ay ~ 1-5. As a result, it is ex- 
pected that the ^^CO emission mainly comes from 
the intermediate part of the density (~ 10^ cm~^) 
and that the C^^O emission comes from the inner 
part of the density 10^ cm~^). 

Figure 1 shows the total-intensity distribution 
of ^^CO (J = 1-0) in the velocity range from Vlsr 
= —30 km s~^ to —10 km s~^, overlaid on an op- 
tical image taken from the Digitized Sky Survey. 
We also show the distribution of OB stars overlaid 
on a ^^CO map and an H^ image in Figure 2. The 
positions of the OB stars are taken from the SIM- 
BAD database. We find ~ 80 0-type stars includ- 
ing five 03 stars, ~400 B stars, and 7 Wolf-Rayet 
stars in the observed area. The H^ image is re- 
produced from the Southern H„ Sky Survey Atlas 
(SHASSA; Gaustad et al. 2001). It is clearly seen 
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that the ovcraU distribution of the ^^CO emission 
shows an anti-correlation with the optical objects 
such as stellar clusters and H ii regions. Most of 
the 0-type stars are located near the rj Car neb- 
ula at (/, b) ~ (287.5, —0.7), creating the huge 
Ha emission nebulosity, while small H ii regions 
are made by a few O stars such as NGC 3293 at 
(285.9, 0.1) and Gum 31 at (286.2, -0.2). On the 
other hand, there is a good correlation of the ^^CO 
emission with the optical obscuration around the 
r] Car nebula; the northern cloud at (287.3, —0.7), 
the southern cloud at (287.7. —0.7), and the giant 
pillar at (288.0, -1.1). 

In the ^^CO emission, a dozen intense peaks 
are distributed within more extended diffuse emis- 
sion. The total molecular mass traced in the ^•^CO 
emission is ^ 3.5 x 10^ Mq, which is estimated 
in the following manner: We assumed an empiri- 
cal relation that the molecular hydrogen column 
density, A^(H2) is proportional to the ^^CO in- 
tegrated intensity, /(^^CO). The conversion fac- 
tor is called the X-factor [X = 7V(H2)//(i2CO)], 
and the value lies in the range (1-3) x 10^° cm~^ 
(K km s~^)~^ toward the clouds in the Galac- 
tic disk, as estimated by the virial mass, 7-ray 
emission, and so on (e.g., Bloemen et al. 1986; 
Solomon et al. 1987; Strong et al. 1988; Bertsch et 
al. 1993). Here we adopted X = 1.6 x 10^° cm'^ 
(K km s~^)~^ as derived from observations with 
the Energetic Gamma-Ray Experiment Telescope 
(EGRET; Hunter et al. 1997). The total mass of 
the ^^CO emission is estimated by using 

M = ixmn 'Yj D'^nN{Il2+B.e) = 2.8mH X] D'^^N{H2), 

(1) 

where D is the distance, Q is the solid angle sub- 
tended by the effective beam size (2' x 2' or 4' 
X 4'), rriH is the proton mass, and is the mean 
molecular weight. The factor of 2.8 is introduced 
by taking into account the contribution of helium 
(one He atom for every five H2 molecules), which is 
the same assumption as /x = 2.33. The summation 
is performed over the observed points within the 3- 
sigma contour level (5 K km s~^) of the integrated 
intensity. The distance of the rj Car GMC is not 
well determined. Grabelsky et al. (1988) used D 
= 2.7 kpc, which is estimated from the distances 
of associated star clusters. Feinstein (1995) sum- 
marized the photometric data of the open clusters 
distributed over the entire ry Car GMC, and ob- 



tained the average value of 2.5 kpc. In the vicin- 
ity of 77 Car, values in the range 2.2-2.8 kpc are 
reported for the distances of Tr 14 and Tr 16 (e.g., 
Davidson & Humphreys 1997). In this paper, we 
adopt the frequently quoted value of 2.5 kpc as the 
distance of the 77 Car GMC (Feinstein 1995). The 
uncertainty in determining the total mass is esti- 
mated to be ~ 70% from the equation (1), since 
the distance D and the conversion factor X in- 
clude ^ 10% and ^ 50% uncertainty, respectively. 
The total mass obtained in this manner is some- 
what smaller than that derived in Grabelsky et 
al. (1988), 6.7 x 10'^ Mq (cloud 11 [287.5-0.5] in 
their catalog), who observed the Carina arm with 
a 1.2-meter telescope at Cerro Tololo. The dif- 
ference in the total mass basically originates from 
the different assumptions on the X-factor and the 
distance used in deriving the mass; they used X = 
2.8 X 10^0 cm-2 (-j^ j^^^ s-i)-i and D = 2.7 kpc. 
When we use the same assumption as in Grabel- 
sky et al. (1988), 7.2 x 10^ Mq is obtained, which 
is in good agreement with the value derived by 
Grabelsky et al. (1988). 

Clumpy structure is more clearly seen in the 
13C0 (J = 1-0) emission (Fig. 3). From the ^^CO 
distribution, we can divide the rj Car GMC into 7 
regions, which are clearly separated at the 3-sigma 
contour level (1.35 K km s~^). The boundaries of 
the regions are similar to those in a previous study 
by Zhang et al. (2001b), which is based on the 
^^CO (J = 1-0) map of 7f5 grid spacing with an 
8f8 beam, but our definition traces more detailed 
structure. Moreover we identified one more region 
(region 7). This is mainly because the spatial res- 
olution in our study is ~ 3 times higher than that 
in Zhang et al. (2001b) and because we used ^^CO 
emission instead of ^^CO. 

The physical properties of the ^"'CO gas arc 
estimated in the following manner assuming lo- 
cal thermodynamical equilibrium (LTE): The ex- 
citation temperature, Tox, was estimated from 
(^^CO) at the peak position of each C^^O core (see 
§ 3.2) by using 



ln{l + 5.53/[T^(i2CO)(K)-h 0.819]}^ 

(2) 

We used the average value of the Tex of each of 
the C^^O cores within the region as the Tex of 
the region. In order to derive the ^^CO column 
density, we divide each spectrum into 0.1 km s~^ 
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bins, calculate the column density within each bin, 
and sum them up within an LSR velocity range 
from —30 km to —10 km s~^. The optical 
depth in each bin, T13 {V), is calculated by the 
following equation: 



■In 1 



5.29{Ji3[Tex(K)]- 0.164} 



(3) 

where T13 (V) is the average temperature of the 
^^CO spectrum in each bin in the unit of Kelvin, 
and Ji3[T(K)] = l/exp[5.29/T(K)] - 1. The ^^CO 
column density is estimated from 



Ni3 = 2.42x10^ 



0.1(kms-i)Ti3(T/)Te,(K) , _2, 



V 



exp[-5.29/Tex(K)] 



-(cm 
(4) 

The ratio A''(H2)/iV(^'''C0) is assumed to be 5 x 
10^ (Dickman 1978). The mass is estimated in the 
same manner as for the ^^CO. The total molecular 
mass traced in the ^^CO emission is ^ 1.3 x 10^ 
Mq, which is 38 % of the ^^CO mass. 

The spatial distribution of the C^^O (J = 1- 
0) integrated intensity is shown in Figure 4. The 
distribution of the C^^O gas is more compact than 
those of the ^^CO and ^^CO gas. The total molec- 
ular mass traced in the C^^O emission is ~ 5.8 x 
10"^ Mq (see § 3.2), which is 44 % of the ^^CO 
mass. 

The mass of each region traced in ^^CO, ^^CO, 
and C^^O is summarized in Table 1. The fraction 
of the ^^CO mass to the ^^CO mass ranges from 
36-68 %, with the largest value of 68 % in region 
1. The fraction of the C^^O mass to the ^-''CO 
mass ranges from 21-72 %, with the largest value 
of 72 % again in regions 1 and 5. 

We note that CO emission with Vlsr ^ — 10 km 
s~^ was detected toward some areas around {I, b) 
= (286.5, —0.2), although these velocity compo- 
nents are not shown in the figures. These velocity 
components come from distant clouds (clouds on 
the far side of the Carina arm ai D = 7.5 kpc, Gra- 
belsky ot al. 1988), not associated with the rj Car 
CMC. We will discuss these clouds in a separate 
paper (Yonekura et al. 2005 in preparation). 

3.2. Identification of C^^O Cores and their 
Physical Properties 

In order to investigate the physical properties 
of the dense regions, we define a C^*0 core in 



the same manner as that adopted by Onishi et 
al. (1996); (1) find a peak-intensity position, (2) 
draw a contour at the half level of the peak inten- 
sity, (3) identify a core unless previously identified 
cores exist within the half- level contour, (4) find 
the next intensity peak outside the core, (5) re- 
peat the procedure after (2) until the peak inten- 
sity falls down below 6-sigma level (= 1.0 K km 

We finally identified 15 C^^O cores. In Figure 5 
we show the boundaries of each core, and their ob- 
served properties are listed in Table 2. The abso- 
lute peak temperature, (C^^O), the radial ve- 
locity of the C^^^O emission, Vlsr, and the FWHM 
line width, AV, are derived from a single Gaussian 
fitting for a line profile at the peak position of the 
core. There is a variation in the ratio of the C^^O 
integrated intensity, /(C^^O), to the ^^CO inte- 
grated intensity, /(^''CO); throe prominent C-'^^O 
cores, whose C^^O intensities at the peak position 
are larger than ~ 2 K km s~^ (nos. 5, 6, and 13), 
possess only moderate ^"^CO intensity among the 
dozens of ^^CO local peaks. This indicates that 
^^CO observation is not necessarily suited to trace 
the high-density regions. 

The physical properties of the C^^O cores were 
estimated using nearly the same procedure as for 
-'^^CO. Differences between the two procedures are 
the following: The optical depth of the C^*0 line 
in each bin, tis{V), is calculated by the following 
equation: 



In 1- 



TisiV) 



5.27{Ji8[Tex(K)] - 0.166} 



(5) 

where Tis{V) is the average temperature of 
the C"'^*0 spectrum in each bin in Kelvin, and 
Ji8[T(K)] = l/exp[5.27/T(K)] - 1. The C^^O 
column density is estimated from 



(6) 

The ratio N{R2)/N{C'^^0) is assumed to be 6 x 
10^ (Prerking, Langer, & Wilson 1982). The mass 
of a C^'^^O core is derived by summing over the 
observed points within the core. The radius of a 
core, r, was calculated from 



(7) 
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where S is the area inside the core. The radius of 
small cores is overestimated by the beam-size ef- 
fect. The relation between the observed radius, r, 
and the deconvolved radius, rdec, is approximated 
by the relation 



beam size 



We also calculated the virial mass, Myir, 



Myj, 

Mr, 



209 — 



pc/ \ km s~^ 



(8) 



(9) 



with AVcomp defined as the FWHM line width 
of the composite profile derived by using a sin- 
gle Gaussian fitting, where the composite profile 
was obtained by averaging all the spectra within 
the core. The average H2 number density of the 
core, n(H2), was derived by dividing M by the vol- 
ume of the core, assuming that the core is spherical 
with radius r. The physical properties of the C^^O 
cores thus obtained are listed in Table 3 along with 
the average and the median values. 

Here we compare the physical properties of the 
cores in the 77 Car GMC to those of massive-star- 
forming cores observed with the NANTEN and 
the 4-m radio telescopes of Nagoya University: 
Orion B (Aoyama et al. 2001), Orion A (Naga- 
hama 1997), Vela C (Yamaguchi et al. 1999a), 
Cepheus 0B3 (Yu, Nagahama, & Fukui 1996), 
S35/S37 (Saito et al. 1999), and Centaurus (Saito 
et al. 2001). All these studies assume LTE in de- 
riving physical properties and use the same as- 
sumption on abundance. Since the definition of a 
C^^O core in Orion A, Vela C, and Cepheus 0B3 is 
different from the one in this work, only values at 
the peak position (i.e., AV and 7V(H2)) are com- 
pared for these regions. In order to be consistent 
for comparison, we used (Mvir/M©) = 209 (r/pc) 
(Ay /km s~^)^ instead of equation (9), and we 
excluded cores with serious beam-size effect, i.e., 
fdec/i' < 0.8, in calculating the average and me- 
dian values of r, n(H2), and M^i^/M according to 
Saito et al. (2001). Hereafter, we call the cores 
with rdec/r > 0.8 as 'resolved' cores. The aver- 
age and median values in each region are listed 
in Table 4. The line width, AV, of the cores in 
the 7] Car GMC ranges from 2.2 to 5.4 km 
with an average value of 3.2 km s~^. This is the 
largest among the sample. The peak column den- 



sity ranges from 0.6-3.5 x 10^^ cm~^ with an av- 
erage of 1.3 X 10^^ cm~^, which is a typical value 
for the GMCs listed here. The core mass ranges 
from 450 to 6100 Mq with an average of 2600 M©. 
It is difficult to compare the absolute values of the 
mass, since there is a dependence on the distance. 
The radius of the resolved cores ranges from 1.6 
to 3.1 pc with an average of 2.2 pc. The radius 
is also dependent on the distance. The average 
number density, n(H2), ranges from 0.5-3.4 x 10'^ 
cm~'^ with an average of 1.2 x 10'^ cm~^, which is 
about one order smaller than the expected value 
for a C^^O core, ^ 10** cm~^. This is mainly due 
to the beam dilution; the resolution of the tele- 
scope is ~2 pc at the distance of the r] Car GMC, 
which is larger than the typical size of the C'^^O 
cores in Orion, ~ 0.5 pc. The ratio of the virial 
mass, Mvir, to the LTE mass, M, ranges from 0.5 
to 2.9 with an average of 1.6, which is the second 
largest value among the whole sample. 

To summarize, the cores in the rj Car GMC 
are characterized by larger values of AV and 
Myir/M on average than those in other massive- 
star-forming regions. 

3.3. Star-Formation Activity in Dense 
Cores: Existence of Large Numbers 
of Non-Star-Forming Cores 

In this subsection, we investigate the star- 
formation activity in C^^O cores in the rj Car 
GMC. The IRAS point source catalog is a use- 
ful database to search for young stellar objects 
(YSOs) embedded in molecular clouds. In order to 
select candidates for recently formed stars among 
the IRAS point sources, we searched for sources 
having cold infrared spectra defined at 25 and 60 
/um, i.e., sources with higher flux densities at 60 
(ITU than at 25 /im with a data quality better than 
or equal to 2 in these two bands. We excluded 
objects from the selection that have already been 
identified as planetary nebulae or external galax- 
ies in the IRAS point source catalog. Among the 
316 IRAS point sources in the observed region, 
73 were found to satisfy the criteria, and only 
two sources are located inside the boundaries of 
C^^O cores, i.e., within the half-intensity contours: 
IRAS 10365-5803 and IRAS 10361-5830 are as- 
sociated with C^^O cores nos. 5 and 6, respectively 
(Figs. 4 and 5). The luminosities of these sources 
are ~30000 Lq (IRAS 10365-5803) and ~21000 
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Lq (IRAS 10361—5830), indicating that these are 
candidates for massive stars or star clusters. The 
properties of these IRAS point sources are summa- 
rized in Table 5. We also searched for YSOs from 
the MSX point source catalog. We searched for 
sources with a data quality better than or equal 
to 3 in Band C 12/im), D 14^m), and E 
(~ 21/xm). Among the 2799 MSX point sources in 
the observed region, 147 were found to satisfy the 
criteria, and 8 sources are associated with 5 C^^O 
cores (nos. 5, 6, 9, 10, and 12; Fig. 6). The prop- 
erties of these MSX point sources are summarized 
in Table 6. 

The detection limit of the IRAS survey toward 
the T] Car CMC is estimated to be ~ 1000 Lq, 
which we derived from the distribution of the lumi- 
nosity of the above 73 YSO-candidate IRAS point 
sources on the assumption that these sources are 
located at the same distance as the ry Car CMC. 
Therefore, we can conclude that massive star for- 
mation has not occured in the cores without IRAS 
point sources, although we can not exclude the ex- 
istence of sources at an evolutionary stage prior 
to the formation of a massive star, which may 
not have enough emission at infrared wavelengths 
(e.g., Sridharan et al. 2002). 

Along with the result in the next subsection, 6 
C^^O cores (nos. 5, 6, 9, 10, 12, and 15) out of 
15 (= 40%) are experienced star formation, and 
at least 2 of 15 (= 13 %) are massive-star forming 
cores in the 77 Car GMC, although the propor- 
tion of star-forming cores may vary when high- 
spatial-resolution molecular line observations and 
more sensitive infrared and (sub)millimeter obser- 
vations are made. It should be noted that as many 
as 13 out of the 15 C^*0 cores, i.e., 87 %, are not 
associated with luminous IRAS point sources [L 
> lO'' Lq). The proportion, 13 %, is much lower 
than in other massive star forming regions at dis- 
tances within a factor of two of the Car GMC: 
28 % in S35/S37 and 86 % in Centaurus III (Saito 
et al. 1999, 2001). We will discuss the inactiveness 
in detail later in § 4.3 

3.4. Search for High-Density Regions and 
Molecular Outflows 

In this subsection, we present the results from 
the search for high-density regions and molecular 
outflows, in order to reveal the physical proper- 
ties and the evolutionary stages of dense cores, es- 



pecially those without IRAS point sources, more 
precisely. 

3.4-1- Identification of H^^CO^ Cores and their 

Physical Properties 

In order to detect high-density regions within 
C^*0 cores, we have observed 14 C^*0 cores (all 
except for C^^O core no. 4) in Hi3C0+ (J = 1- 
0) emission, which traces a typical density of ~ 
10^ cm~^(e.g., Mizuno et al. 1994; Onishi et al. 
2002). Therefore, the C^^O cores without IRAS 
point sources, which are detected in H^^CO+ emis- 
sion, can be considered as good candidates for the 
sites of future massive star formation. Significant 
H^^CO"'" emission above the 3-sigma noise level 
was detected toward 4 C^^O cores (nos. 5, 6, 9, 
and 15), two of which (nos. 5 and 6) have associ- 
ated IRAS point sources and three of which (nos. 
5, 6, and 9) with MSX point sources (see § 3.3). 
Distributions of H-'^^CO"'" emission are presented 
in Figure 7. In order to study the physical prop- 
erties, we define the H^^CO+ cores in the same 
manner as C^^O cores. As a result, we identified 
3 cores. 

Physical properties of the H^^CO+ cores were 

estimated using the same method as for C^*^0: The 
H^^CO"*" optical depth in each bin, tv, is calcu- 
lated from 



T(H^3C0+)y = -ln( 1 



T(Hi3C0+); 



4.16 [J(Tex) - 0.272] 

(10) 

where J(T) = 1 / [exp(4.16/r) - 1]. We tenta- 
tively estimate Tox of each H^^CO+ core from the 
^^CO (J = 2-1) observations (see § 3.4.2), since 
the beam sizes of both observations are similar. 
The H13C0+ column density, N{n^^CO+), is es- 
timated by using the following equation: 



N{^^^CO+) = 1.92x10^^^ 



0.1(km s-^)T(Hi3C0+)yTex(K) 



l-exp[-4.16/rex(K)] 
(11) 

The II2 column density, A/'(H2), was derived from 
iV(H^''CO+) by making the same assumption as in 
Aoyama et al. (2001); [HC0+]/[H2] = 4.0 x lO"'^ 
and i^c/i^c = 89. The assumption of Tex may 
not be correct, since the density is perhaps smaller 
than the critical density oiB^^CO^ (J = 1-0), i.e., 
J = 1-0 line of H^^GO+ is perhaps sub-thermally 
excited. However, the uncertainty in T^x would 
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not seriously affect the estimate of the column 
density A^: Even if we changed the value of Tex 
by a factor of 2, the estimated value of N changes 
within a factor of 1.9, which is smaller than the un- 
certainty in the HCO"*" abundance (e.g., Bergin 
et al. 1997; J0rgensen, Schoier, & van Dishoeck 
2004). 

The physical properties of the H^'^CO+ cores 
are listed in Table 7. The total mass of the molec- 
ular gas traced by H^'^CO+ in the observed area 
above the 3-sigma level is estimated to be ~ 4300 
Mq, corresponding to ~ 7 % of the C^^O mass. 

It is to be noted that a significant difference is 

found in the behavior of the line width between 
cores with and without IRAS point sources: The 
line widths in the II^^CO+ cores are much smaller 
than those in the C-'^^O cores for the cores without 
IRAS point sources (nos. 9 and 15), whereas the 
line widths in the H^^CO"*" cores are nearly the 
same as or somewhat larger than those in C'^^O 
cores for cores with IRAS point sources (nos. 5 
and 6). The present result can be compared to 
the multitracer, single-cloud line width-size rela- 
tion, which is classified as Type 3 in Goodman et 
al. (1998). The relation AV oc R" with 0.2 < a 
< 0.7 is reported with the tendency that massive 
star-forming regions tend to have lower values of a 
(e.g., Caselli & Myers 1995). The present result is 
consisitent with the above tendency in a sense that 
the value a in massive star-forming cores tends to 
be smaller than in the others, although the trend 
is more extreme: i.e., a < 0. The trend might 
be indicating that line-width enhancement due to 
the input of turbulence by the activity of star- 
formation, such as jots and molecular outflows, is 
effective especially near the center of the massive 
star-forming cores. However, we have to wait un- 
til statistically more accurate studies can be made 
based on a large sample. 

This section can be summarized as follows: We 
detected H'^'^CO"'" emission from 4 C^^O cores, 
which indicates that high-density gas is presented 
in these cores and that these cores are capable of 
star formation. Indeed, two (nos. 5 and 6) have al- 
ready experienced massive star formation and an- 
other one (no. 9) has experienced less-massive star 
formation. The other C'^^O core (no. 15) may have 
not yet experienced star formation, which makes 
the core the most probable site for future star for- 
mation, though we cannot rule out the possibil- 



ity that less-massive stars are forming inside (see 
§ 3.4.2). 

3.4-2. Molecular Outflows 

It is well known that molecular outflows are 
ubiquitous in low-mass star formation (Lada 1985; 
Fukui et al. 1986; Fukui 1989; Fukui et al. 1993 
and references therein), and they are believed to 
be related to the accretion process. It is also sug- 
gested by these studies that molecular outflows 
arc also common in massive star forming regions. 
Recent systematic surveys for outflows toward re- 
gions of massive-star formation showed a detection 
rate of molecular outflow phenomena as high as 
~90 %, confirming that outflows are also common 
in massive-star forming regions (e.g.. Shepherd & 
ChurchweU 1996: Zhang et al. 2001a; Beuthcr et 
al. 2002b). The formation process of massive stars 
is, however, still controversial; it is not yet known 
whether massive stars are formed by direct accre- 
tion on a central protostar through an accretion 
disk (e.g., McKee & Tan 2003), or by the coa- 
lescence of low- to intermediate-mass protostars 
(Bonnell, Bate, & Zinnecker 1998). At any rate, 
the detection of a molecular outflow is a direct 
indication of on-going star formation. The age 
and the luminosity of the driving source can be 
estimated from the observations, and moreover, 
the kinematic information, such as the interaction 
with parent molecular clouds as well as the sys- 
temic velocity of the driving source, can be used to 
solve the distance ambiguity in cases where more 
than one molecular cloud is located in the same 
line of sight. 

We searched for molecular outflows toward 7 
C^^O cores (nos. 5-9, 14, and 15). We detected 
one outflow as well as three candidates for out- 
flows. Figure 8 displays the '^^CO spectra around 
C^^O core no. 15. In the flgure, both the blue- 
shifted (-35 km s"^ < Vlsb. < -25 km s'^) and 
red-shifted wings (—15 km s"'^ < Vlsr < — 5 km 
s~^) are clearly seen in some spectra, where the ve- 
locity ranges for wing components are determined 
at a 1-(T noise level in the spectra smoothed to 
the velocity resolution of 0.5 km s~^ (T^ ~ 0.8 
K): The blue-shifted component is prominent at 
the offset (Aa, AS) = (-1', 0'), while the red is 
most intense at (—Of 5, 0'). The distribution of 
the ^^CO line intensities integrated between Vlsr 
= —35 and —25 km s~^ (blue component) and be- 
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twccn —15 and —5 (red component) is shown in 
Figure 9. 

Figure 10 shows ^^CO profiles of three outflow 
candidates. Although their spectra show asymme- 
try, or have broad wing features, we cannot make 
a contour map of the wing components, due to the 
existence of multiple velocity components and/or 
contamination from the emission at OFF position. 

Physical parameters of the detected outflow as- 
sociated with the C^^O core no. 15 were estimated 
as follows: The excitation temperature of the high- 
velocity gas was assumed to be the same as that 
of the quiescent gas, and estimated to be 30 K 
from the ^^CO (J = 2-1) observations by using 
the following equation 



ln{l + 11. l/[T*(i2C0, 2-l)(K) + 0.187]}^ 

(12) 

The column densities of ^^CO (J = 2-1) are cal- 
culated under the assumption of LTE using the 
following equation: 

N{CO) = 1.04x10^3 Tex exp (^^^ ^ J T^C^CO)dV, 

(13) 

where the integration is made over the wing com- 
ponent; P is defined from the ^^CO optical depth 
as /3 = [1 — exp(— r)] / r. We assumed r = 2 for 
both the red and blue components, which is a typi- 
cal value for outflows observed in the ^^CO (J = 2- 
1) line (Levreault 1988). The column density of H2 
molecules was estimated using an [H2]/[C0] ratio 
of 1 X lO'' (Frerking, Langer, & Wilson 1982). The 
mass of each component is derived using equation 
(1) with an effective beam size of 30" x 30". The 
summation is performed over the observed points 
within the 3-sigma contour level (5 K km s~^) of 
the intensity integrated over the wing component. 
Sizes of the lobes are defined at the 3-sigma con- 
tour level of the integrated intensity map of the 
outflow. We tentatively defined the position of 
the driving source of this outflow as the center 
of the peak position of each lobe, because no driv- 
ing sources, such as IRAS and MSX point sources, 
have been identified around the outflow yet. The 
radius is defined as the maximum separation be- 
tween the position of the driving source and the 
3-sigma contour of the wing emission. 

Using the mass and the radius of the outflow, 
we estimate the kinetic energy, mechanical lumi- 



nosity, momentum, and dynamical timescale of the 
outflow in the same manner as in Yonekura et al. 
(1998). We calculated physical properties of the 
outflow for the velocity ranges from —35 to —25 
km s~^ for the blue-shifted wing and from —15 to 
—5 km s~^ for the red-shifted wing. 

These values are summarized in Table 8. The 
nature of the driving source is not clear because 
no sources have been identified yet around the out- 
flow. We estimate the upper limit to the luminos- 
ity from the IRAS sky survey atlas to be 3900 
Lq. Another estimate for the luminosity can be 
made by using the empirical relation between the 
mechanical luminosity of the outflow and the bolo- 
metric luminosity of the driving source (Cabrit & 
Bertout 1992); 15 Lq is obtained in this manner. 
There are two possibilities for the nature of the 
driving source; [i] a less-massive star with a lu- 
minosity below the IRAS and/or MSX detection 
limit or [ii] a source at an evolutionary stage prior 
to the formation of a star without enough emission 
in the infrared (e.g., Sridharan et al. 2002). We 
cannot rule out either of the above two possibili- 
ties at the moment, but the former might be more 
probable because the dynamical timescale (~ 3 x 
10** yrs), which will probably be overestimated be- 
cause of the overestimation of size due to low an- 
gular resolution, is relatively long. High-resolution 
CO observations as well as sub-millimeter dust 
continuum observations will reveal the nature of 
the driving source. 

4. DISCUSSION 

4.1. Physical Properties of C^'^O Cores 
and their Relation to Star Formation 

We discuss here the relationship between the 
star-formation activity and the physical properties 
of C^*0 cores, which has become possible thanks 
to the detection of a large number of cores without 
any signs of star formation. 

In the study of C^^O cores in Orion B, Aoyama 
et al. (2001) found that star formation occurs pref- 
erentially in cores with larger AV, N, and M, and 
that the luminosity of the forming stars increases 
with TV. In Centaurus, the luminosity of the form- 
ing stars increases with N and M as well as AV, 
and decreases with Mvir/M (Saito et al. 2001). 
Saito et al. (2001) also found that the luminosity 
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increases with AV by compiling the data in Orion 
A, B, S35, Vela C, and Centaurus, and proposed 
a possibility that the mass of forming stars is de- 
termined by AV , since high mass-accretion rates 
can be achieved in cores with large Ay. However, 
it is controversial whether large AV is a cause or 
effect of massive star formation. 

We divided the C^^O cores into 3 groups: (A) 
cores with IRAS point sources (nos. 5 and 6), (B) 
other star-forming cores (nos. 9, 10, 12 and 15), 
and (C) cores without any signs of star forma- 
tion. The cores in group A are thought to be sites 
of massive star formation (see Table 5 for the lu- 
minosities of the associated IRAS point sources). 
There are two possibilities for the nature of the 
cores in group B; these are either [i] less- massive 
stars with luminosities below the IRAS detection 
limit or [ii] sources at an evolutionary stage prior 
to the formation of stars without enough emission 
in the infrared (see § 3.4.2). Here we tentatively 
treat them as less-massive star- forming cores, fol- 
lowing the argument in § 3.4.2. The cores in group 
C seem to be sites of less-active star formation 
than those in groups A and B, which is inferred 
from the absence of the high-density regions (see 
§ 3.4.1). Table 9 gives the average physical prop- 
erties of these 3 groups as well as the average of all 
cores in the 77 Car GMC. We also calculated the 
average properties of all the star-forming cores, 
i.e., cores in groups A and B. Figure 11 shows a 
series of histograms of the core parameters, [a] line 
width AV^omp, [b] H2 column density A''(H2), [c] 
LTE mass M, [d] radius r, [e] H2 number density 
n(H2), and [f] ratio Myi^/M. Note that only re- 
solved cores are used to make the histograms of r, 
n(H2), and Mvir/M. 

It is clearly seen that iV(H2), M, and n(H2) are 
larger and the ratio Mvir/M is smaller in group A 
than in groups B and C, while the tendency is not 
clear in AT^omp and r (trend 1). Between star- 
forming cores (groups A-l-B) and non-star-forming 
cores (group C), A^(H2) and n(H2) are larger in 
groups A+B than in group C (trend 2). How- 
ever, no significant difference between group B and 
group C is seen (trend 3). These trends (trends 1- 
3) might be indicating the differences in the physi- 
cal properties among the cores as a function of the 
masses of the forming stars. 

Here we compare the present results to those in 
Centaurus (Saito et al. 2001). In Centaurus, the 



luminosity of the forming stars increases with 
and M as well as AV, and decreases with Myir/M. 
This trend is same as the trend 1 in 77 Car GMC 
except AV. Although we have obtained 127 C'^^O 
cores in massive star forming regions so far by us- 
ing the same telescope including 15 in 77 Car GMC, 
we can not make direct comparison at present, 
because of the differences in the detection limit 
and/or the definition of C^^O cores. The results 
of detailed comparion will be published in a sepa- 
rate paper. 

We can conclude that massive star formation 
occurs preferentially in cores with larger N, M, n, 
and smaller Mvir/M in rj Car GMC. High spatial- 
resolution and high sensitivity observations in in- 
frared (e.g., Spitzer, ASTRO-F) as well as the ob- 
servations in high density tracer such as CS (J 
= 7-6), HCO+ (J = 4-3), HCN (J = 4-3) (e.g., 
ASTE), will help to reveal the properties of C^^O 
cores. 

4.2. Comparison with Sub-mm Observa- 
tions 

It is of interest to compare the distribution of 
CO and C^^O (J = 1-0) emissions to that of CO 
(J = 4 3) emission obtained with the Antarctic 
Submillimeter Telescope and Remote Observatory 
(AST/RO; Zhang et al. 2001b), since different ex- 
citation conditions arc required between CO (J 
= 4-3) and CO (J = 1-0). The J = 4 level of 
the CO molecule is 55 K above the ground level, 
and the critical density of the J = 4 3 transition is 

8 X 10** cm~"^ for the optically thin case, whereas 
the J = 1 level of the CO molecule is 5.5 K above 
the ground level, and the critical density of the J — 
1-0 transition is ~ 10^ cm~^. At kinetic tempera- 
ture of 40 K, the J = 4-3 emission of CO is most 
intense at density of ^ 10"'' cm~'^ while the J = 
1-0 intensity is peaked at around 10"^ cm~^ (Gold- 
smith and Langer 1978). Observationally, the den- 
sity effectively emitting these emissions may cor- 
respond to somewhat lower than these values (~ 
10^ cm~^ for the J = 1-0 emission), which can 
be understood if we take into account the optical 
depth effect. Thus, the CO (J = 4-3) line requires 
higher temperature and density than the CO (J = 
1-0) line. However, the overall distribution of the 
molecular gas traced in the CO (J = 1-0) line 
and that of the CO (J = 4-3) line is quite similar 
except the following two regions: 
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(1) C^^O core no. 10 (adjacent to Tr 14) 

The integrated intensities of CO (J = 4-3) and 
CO (J = 1-0) are the strongest in the observed 
region, although the integrated intensity of C^^O 
(J = 1-0) shows only moderate values. A possible 
explanation for this is a high excitation tempera- 
ture due to the star cluster Tr 14: Because of the 
high excitation temperature, the optical depth of 
C^*0 (J = 1-0) remains very small in spite of its 
large column density, resulting in very low inten- 
sity in C^^O (J = 1-0) emission. The NANTEN 
data indicate the excitation temperature is ~ 25 
K, and this value seems to be a lower limit, since 
higher resolution data obtained with a 24" beam 
show the excitation temperature to be as high as 
- 45 K (Brooks 2000). 

(2) C^^O core no. 5 

The intensity of CO (J = 4-3) is very weak, 
whereas CO and C^^O (J = 1-0) emission is quite 
strong. In general, two explanations are possible; 
(1) the temperature is not high enough to excite 
the J = 4 3 emission and/or (2) the density is 
lower than the typical value traced in CO (J = 4- 
3). However, the excitation temperature derived 
from the CO (J = 1-0) observations at NANTEN 
is 28 K, the highest value among the C^^O cores 
in the observed region. The number density, 2800 
cm""^, is also the highest. Moreover, the detec- 
tion of H^'^CO"'" (J = 1-0) emission toward the 
core indicates that a high-density region of > 10^ 
cm~^ exists within the core. Thus, the above two 
hypotheses fail to explain the properties of the 
core. High-resolution observations of both low- J 
and mid- J CO lines may unveil the properties of 
the core. 

4.3. Characteristics of the Star Formation 

in the ry Car CMC 

In this subsection, we will describe the charac- 
teristics of the star formation in the rj Car GMC. 
The proportion of massive-star forming cores, i.e., 
cores with IRAS point sources with L > 10'^ Lq, 
13 % (= 2/15), is much lower than in massive star 
forming regions at distances within a factor of two 
of the T] Car GMC; 28 % in S35/S37 and 86 % 
in Centaurus HI (Saito et al. 1999, 2001). This is 
consistent with the fact that the ratio of M^i^/M 
is large (see § 3.2), i.e., most of the cores are not 
gravitationally bound. The large line widths of 
the cores (~ 3 km s~^) compared to the other 



GMCs 2 km s"'^, sec § 3.2) also support this. 
A large amount of turbulence in the r] Car GMC 
may prevent star formation at the present time. 

Let us then consider the origin of the large tur- 
bulence in the ry Car GMC. The kinetic energy due 
to the turbulence is roughly estimated from _Ekin 
= i M V^^^^. By using the ^^CO mass M = 3.5 x 
10^ Mq, and the average line width of the C^^O 
cores Vturb = 3 km s"\ Skin ~ 3.4 x lO"'^ ergs 
is obtained as the kinematic energy of the q Car 
GMC. The momentum P = M Vturb is estimated 
to be 1.0 X 10^ Mq km s~^. In the absence of any 
mechanisms to prevent it, the turbulence would 
dissipate on the crossing time-scale of the cloud, 
tcr ~ ^?cioud/AKioud (sec Vazqucz-Semadeni et al. 
2000 for a review), which corresponds to ~ 15 Myr 
in the case of the r] Car GMC. Thus the presence 
of a large amount of turbulence means either that 
the turbulence is somewhat constantly injected by 
some mechanisms, or that a large amount of tur- 
bulence is still retained although the turbulence 
is dissipating because the GMC is young and/or 
the initial amount of the turbulence is large. In 
the foUowings, we examine four possibilities: (1) 
Turbulence injection by mass outflows from young 
stellar objects, (2) by supernova remnants, (3) by 
stellar winds from massive stars, and (4) turbu- 
lence is pre-existing. 

(1) Turbulence injection by mass outflows from 
young stellar objects: The kinetic energy and the 
momentum of such outflows range from E = 10^'^- 
10^'^ ergs and from P = 0.1 1000 Mq km s"!, 
respectively (Fukui et al. 1993). Thus, even if we 
assume that all outflows have E — 10''^ ergs and 
P = 1000 Mq km s^^ and that all the energy 
and momentum of such outflows arc converted into 
turbulence (i.e., conversion coefficient e = 1) at 
least 1000 YSOs are needed to supply the present 
turbulence of the rj Car GMC. 

Here we estimate the contribution of outflows 

from low-mass YSOs, since most of the stars are 
formed as clusters in massive star forming regions. 
Observationally, it is known that the kinetic en- 
ergy and the momentum of outflows are corre- 
lated with the bolometric luminosity of the driving 
source, Lboi (Wu et al. 2004). Least-squares fits to 
their data yield P = lO^O ''*^ LjJ'^f and E = lO^i *^^ 
-^boi*- If assume that the relation between the 
mass and the luminosity of YSOs as L = M" with 
a = 3.45 (Allen 1972) and that the initial mass 
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function (IMF) of YSOs as Scale's IMF (dn/dM 
= M-^ b = 2.7, Scalo 1986) with the cutoff at 
1 Mq, 1.6 X 10^ YSOs are needed to supply the 
present turbulence of the 77 Car GMC. The number 
of YSOs largely depends on the cutoff mass of the 
IMF and the indices a and b, however, 3 x 10^ are 
still needed for the IMF with the cutoff at 3 Mq 
with a = 4 and 6 = 2. Although these estimate 
are quite rough, the number required to supply the 
present turbulence of the rj Car GMC is extremely 
large compared to the number of YSOs formed in 
a single GMC in any assumption. Therefore, we 
can conclude that the turbulence injection by mass 
outflows from YSOs are not important. 

(2) Turbulence injection by supernova rem- 
nants: The kinetic energy and the initial speed 
of ejecta of a single supernova remnant are typi- 
cally iJkin = 10^^ ergs and V = 10^ km s"^. If 
we assume the mass of the ejecta to be M = 3 
Mq, P = 3 X 10^ Mq km s'^ can be provided 
by a single supernova. Thus, even if we assume e 
^ 1, at least ~30 supernovae are needed to sup- 
ply the present turbulence of the rj Car GMC, 
which is too large when we consider the ages of 
the associated clusters, 3-7 Myr (e.g., Feinstein 
1995). The fact that no supernova remnants are 
known to exist (Whiteoak 1994) supports this, al- 
though possible non-thermal radio emission is re- 
ported (Jones 1973; Tateyama, Strauss, & Kauf- 
mann 1991). Fukui et al. (1999) suggested the 
possible influence of the Carina flare supershell, 
which may have been formed by more than 20 
supernovae located 100 pc away from the rj Car 
GMC. The momentum supplied to the 77 Car GMC 
from the Carina flare supershell through the cross 
section of the r? Car GMC is estimated in the fol- 
lowing manner: The total momentmn of the Ca- 
rina flare supershell is assumed to be 20 times the 
typical momentum of a single supernova. We also 
assume the depth of the GMC to be the same as 
the projected extent of the GMC, ~ 130 pc. The 
momentum supplied through the cross section of 
the rj Car GMC is estimated to be P = 3 x 10'' x 
20 X (130 X 130)/(47r x 100^) = 8 x lO'' Mq km 
s~^ on the assumption of £ = 1, and thus it fails 
to explain the origin of the turbulence. 

(3) Turbulence injection by stellar winds from 
massive stars: The terminal velocity of the stellar 
wind and the mass-loss rate from 77 Car are V = 
700 km s~i and M = 5 x lO""* Mq yr-\ respec- 



tively (van Boekel et al. 2003; Smith et al. 2003). 
During the time scale of 10^ yr, 77 Car will be able 
to provide Ekin = 2 x 10^^ ergs and P = 3 x 10^ 
Mq km s^^, which is comparable to the energy 
and momentum needed if the conversion efficiency 
is ~ 100 %. However, the effect of the stellar wind 
will be confined to a small region around rj Car; 
a radius of ~ 25 pc is derived from the models of 
Castor, McCray, & Weaver (1975) and Weaver et 
al. (1977) by taking the typical density of a molec- 
ular cloud, 100 cm~'^, and thus injection from 77 
Car alone is not sufficient to explain the observed 
turbulence. 

There are 7 WR stars and ^80 0-type stars 
around the 77 Car GMC (see Fig. 2). A typical 
value of the terminal velocity of the stellar wind 
and the mass-loss rate of a WR star is ~ 2000 
km s"-' and M = 4 x 10~^ Mq yr~^ (van der 
Hucht 2001). The total energy and momentum 
from these WR stars in 10^ yr are estimated to 
be i^kin = 1.1 X 10^2 gi-gs and P = 6 X 10^ Mq 
km s~^, respectively. For an 0-type star, typical 
values are V - 1500 -2500 km s'^ and M ~ 10"^- 
10""^ Mq yr^i (e.g., Markova et al. 2004). Over a 
time scale of 10^ yr, one 0-type star will be able 
to provide iJkm ~ 2.3 x 10''*-6.3 x 10^° ergs and 
P ~ 1.5 X 10^-2.5 X 10^ Mq km s'^ Thus about 
one hundred O stars can provide enough E'kin and 
P on the assumption of £ = 1. It may therefore 
be possible that somewhat uniformly distributed 
WR stars and 0-type stars account for the large 
turbulence, although the estimate of e is highly 
uncertain. 

(4) Turbulence is pre-existing: Bertoldi & Mc- 
Kee (1992) argued that massive cores are magneti- 
cally supercritical. This conclusion is extended to 
low-mass cores by Nakano (1998). Magnetically 
supercritical cores are supported by turbulence 
against the pressure of the surrounding medium 
Pg. When Pg exceeds a critical value Per, the core 
cannot be in magnetohydrostatic equilibrium and 
collapses. Nakano (1998) suggested that dissipa- 
tion of turbulence is the most important process 
in reducing Per. The condition Per < Pg can, how- 
ever, be achieved without dissipating turbulence 
if Pg is large. Thus, star formation in molecular 
clouds with large turbulence is possible. Here we 
investigate the possibility that the 77 Car GMC still 
retains a large amount of turbulence although that 
turbulence is dissipating, i.e., the 77 Car GMC is 
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not a relaxed system. There are two possibilities; 
the 1] Car GMC is young and/or it has an initially 
large turbulence. The former may be inconsistent 
with the fact that the ages of the associated clus- 
ters are not extremely small, 3-7 Myr (e.g., Fe- 
instein 1995), and thus we concentrate the latter 
possibility. 

Concerning molecular cloud formation, various 
theories have been proposed so far; cloud forma- 
tion by spontaneous instabilities, by the compres- 
sion of diffuse interstellar medium, and by random 
coalescence of existing clouds (see Elmcgrcen 1993 
for a review). The most realistic scenario may be 
a combination of the above processes. The origin 
of the turbulence in molecular clouds is not yet 
understood, but possible formation mechanisms 
of a molecular cloud with large turbulence have 
been proposed: For example, Koyama & Inutsuka 
(2000) investigated the propagation of a shock 
wave into the interstellar medium and suggested 
that the postshock region collapses into a cold 
layer through thermal instability. Their recent 
two-dimensional calculation shows that the layer 
breaks up into small cloudlets, in which a large 
amount of molecules are formed, in several Myr 
and that the velocity dispersion of the cloudlets 
is typically several km s~^ (Koyama & Inutsuka 
2002). The fragments then coalesce and become 
larger clouds with supersonic velocity dispersions. 
In the case of the rj Car GMC, multiple super- 
nova explosions may have contributed to a certain 
degree; Fukui et al. (1999) suggested that ~20 su- 
pernova explosions may have occurred ~ 100 pc 
away from the rj Car GMC over the last 2 x lO'' 
yr. These supernova explosions may have trig- 
gered the formation of the r] Car GMC as well as 
the formation of massive stars within the GMC. 

Although the above estimate is quite rough, we 
suggest the possibility that a large amount of tur- 
bulence was supplied when the GMC was formed, 
and is now dissipating. To identify the forma- 
tion mechanism for the rj Car GMC is beyond the 
scope of the paper, but compression due to super- 
nova explosions in the Carina flare supershell may 
have played an important role. Turbulent injec- 
tion from stellar winds, supernova remnants, and 
outflows originating in the stars formed within the 
GMC may also contribute, supplying turbulence 
to some degree. 



5. SUMMARY 

We carried out an unbiased survey for massive 
dense cores in the giant molecular cloud associ- 
ated with T] Carinae with the NANTEN telescope 
in i^CO, and C^^O J = 1-0 emission hues. 

The purpose of the present study is to obtain a 
sample of massive-star/ess-corc candidates which 
are possible sites for future massive star formation. 
With moderate spatial resolution, 2'.7, our obser- 
vations covered the entire r] Car GMC (3x2 deg^) 
including regions farther from rj Car, where no pre- 
vious observations had been made except coarse 
large-scale surveys. We also made high-resolution 
observations in Hi3C0+ (J = 1-0) with the MO- 
PRA 22-m millimeter telescope and searched for 
molecular outflows in ^^CO (J = 2-1) toward se- 
lected regions with the ASTE 10-m sub- millimeter 
telescope. The main results are summarized as fol- 
lows: 

(1) We identified 15 C^^O cores, whose typi- 
cal line width AV^ompj radius r, mass M, column 
density Af(H2), and average number density n(H.2) 
were 3.3 km s~ \ 2.2 pc, 2.6x10^ Mq, 1.3x10^2 
cm~^, and 1.2x10^ cm~"^, respectively. 

(2) Two of the 15 cores are associated with 
IRAS point sources whose luminosities are larger 
than 10'* Lq, which indicates that massive star 
formation is occuring within these cores. Five 
cores including the two with IRAS sources are as- 
sociated with MSX point sources. We detected 
Hi''CO+ (J = 1-0) emission toward 4 C^®0 cores, 
two of which are associated with IRAS and MSX 
point sources, another one is associated only with 
an MSX point source, and the other is associated 
with neither IRAS nor MSX point sources. The 
core with neither IRAS nor MSX point sources 
shows the presence of a bipolar molecular outflow 
in ^^CO (J = 2-1), which indicates that star for- 
mation is also occming in the core, and the other 
three of the four H^'^CO+ detections show wing- 
like emission. In total, six C^^O cores out of 15 (= 
40%) are experienced star formation, and at least 
2 of 15 (= 13 %) are massive-star forming cores in 
the r] Car GMC. 

(3) We found that massive star formation oc- 
curs preferentially in cores with larger A''(H2), M, 
n(H2), and smaller ratio of M^-„/M. We also 
found that the cores in the 77 Car GMC are charac- 
terized by large AF and Mvir/ M on average com- 
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pared to the cores in other GMCs observed with 
the same telescope. These properties of the cores 
may account for the fact that as much as 60-87 % 
of the cores do not show any signs of massive star 
formation. 

(4) We investigated the origin of a large amount 
of turbulence in the rj Car GMC. We found that 
turbulence injection from stellar winds, molecu- 
lar outflows, and supernova remnants which origi- 
nated from stars formed within the GMC, are not 
enough to explain the existing turbulence. We 
propose the possibility that the large turbulence 
was pre-existing when the GMC was formed, and 
is now dissipating. Mechanisms such as multiple 
supernova explosions in the Carina flare supershell 
may have contributed to form a GMC with a large 
amount of turbulence. 
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Fig. 1. — Integrated intensity map of the ^^CO (J = 1-0) emission in the velocity range VhSR = —30 km 
s~^ to —10 km s~^, overlaid on the optical image taken from the Digitized Sky Survey. The contour levels 
are every 20 K km s~^, starting from 5 K km s~^ (3 sigma). 
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Fig. 2. — Distribution of OB stars is compared with the distribution of ^^CO gas, overlaid on the Hq image 
reproduced from SHASSA. YeUow, green, red, and blue crosses denote 77 Car, WR stars, O stars, and B 
stars, respectively. The sizes of the crosses denote the luminosity of the stars. The contour levels are the 
same as in Figure 1. 
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Fig. 3. — Total intensity map of -'^'^CO (J = 1-0) integrated over the velocity range Flsr = —30 km s ""^ to 
— 10 km s~^. The contour levels are every 1.35 K km s~^ (3 sigma), starting from 1.35 K km s~^. 
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Fig. 4. — Total intensity map of C^^O (J = 1-0) integrated over the velocity range Vlsr = —30 km to 

— 10 km s~^. The contour levels are every 0.5 K km (3 sigma), starting from 0.5 K km s~^. IRAS point 
sources selected as candidates for protostars are also shown as crosses. The sizes of the signs are proportional 
to the luminosity. Two IRAS point sources which are associated with C^^O cores are labelled. 
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Fig. 5. — Positions of the cores indicated on a C^*0 integrated intensity map. The contour levels are the 
same as in Figure 4. Bold contours represent the boundaries of each core. 
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Fig. 6. — MSX point sources selected as candidates for protostars are shown as crosses. Eight MSX point 
sources which are associated with C^*0 cores are labehed. The contour levels are the same as in Figure 4. 
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Fig. 7. — (a) Distribution of the H^^CO+ (J = 1-0) emission toward C^^O core no. 5. The bold contours 

represent H^'^CO^ intensity integrated over the velocity range —23 km s^^ < T^lsr < —16 km s^^. Contour 
levels are every 0.4 K km s~^ (3 sigma), starting from 0.4 K km s~^. The thin contours represent the 
distribution of the C^^O core observed by NANTEN, with contour levels of 50 %, 60 %, 70 %, 80 %, and 90 
% of the peak intensity. The coordinates arc shown in offsets from the reference position (a2000 = lO'^ 38™ 
41!8, (52000 = —58° 20' 4"). The filled circle and the cross denote IRAS and MSX point sources, respectively, 
(b) Same as Fig. 7a but for C^^O core no. 6. The velocity range, the 3-sigma level, and the reference position 
are -30 km s^^ < Vlsr < -20 km s^^, 0.42 K km s"i, and {10^ 37™ 58!0, -58° 46' 40"), respectively, (c) 
Same as Fig. 7a but for C^^O core no. 9. The velocity range, the 3-sigma level, and the reference position 
are -20 km s'^ < Vlsr < -16 km s'S 0.35 K km s-\ and (lO"^ 42™ 53!0, -59° 25' 44"), respectively, (d) 
Same as Fig. 7a but for C^^O core no. 15. The velocity range, the 3-sigma level, and the reference position 
are -20 km s"! < Vlsr < -16 km s'^, 0.45 K km s-\ and (10^ 47™ 57;5, -60° 26' 26"), respectively. 
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Fig. 8. — ^^CO spectra observed toward C^^O core no. 15. The offsets from the position (a2000) '^2000) = 
(10^ 47™ 57^5, —60° 26' 26") are shown in parentheses in the unit of arcminutes. The spectra are smoothed 
to the velocity resolution of 0.5 km s~^. 
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Fig. 9. — Distribution of the molecular outflow within the C^^O core no. 15. The solid contours represent 
the ^^CO intensity of the blue component integrated between Vlsr = —35 and —25 km s^^, and the dashed 
contours show that of the red component between Vlsr = —15 and —5 km s^^. Contours are every 5 K km 
s^^ (3 sigma) from 5 K km s^^ as the lowest. 
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Fig. 10. — -^^CO profiles for the outflow candidates. The positions are (a) a20oo = lO'^ 38™ 34!2, i52ooo = 
-58° 19' 34" {{Aa, A6) = {-!', 0.'5) offset from the reference position in Fig. 7a), (b) a20oo = 10^ 37"" 54!1, 
(^2000 = -58° 46' 40" ((Aa, Ad) = (-0f5, 0') offset from the reference position in Fig. 76), and (c) a20oo = 
IQh 42^ 49!0, (^2000 = -59° 25' 44" ((Aa, AS) = {-O'.b, 0') offset from the reference position in Fig. 7c). 
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Fig. 11. — Histograms of (a) line width, (b) H2 column density, (c) mass, (d) radius, (e) average H2 number 
density, and (f) ratio M^-^jM of cores in the 3 groups. The dark-shaded and light-shaded regions indicate 
massive-star forming cores (group A) and less- massive-star forming cores (group B), respectively. The open 
regions indicate less- active star- forming cores (group C). Note that only resolved cores are used to make the 
histograms of r, n(H2), and M^-^jM. 
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Table 1 

Comparison between ^^CO, ^^CO, and C^*0 Masses in each Region 



Region 




Imax 


6mi„ 


bmax 


Ml2* 


Mi3» 


Mi8» 


M13/M12 


M18/M13 




C^*0 Cores 


Associated Objects 






n 


n 


n 












(K) 






1 


287.800 


288.400 


-1.467 


-0.933 


2.2 


1.5 


1.07 


0.68 


0.72 


20 


13, 15 


giant pillar 


2 


287.600 


287.900 


-0.900 


-0.400 


1.7 


0.8 


0.28 


0.44 


0.37 


20 


12 


r/ Car, southern cloud 


3 


286.867 


287.567 


-1.000 


0.067 


9.8 


3.8 


0.80 


0.39 


0.21 


23 


8, 9, 10, 11 


rj Car, northern cloud 


4 


286.033 


286.967 


-0.800 


-0.033 


8.0 


3.0 


1.00 


0.38 


0.33 


21 


6, 7 


Gum 31 


5 


285.800 


286.533 


-0.200 


0.333 


6.6 


3.2 


2.29 


0.48 


0.72 


21 


2, 3, 4, 5 


Gum 31, NGC 3293 


6 


285.700 


286.000 


-0.800 


-0.300 


2.0 


0.7 


0.27 


0.36 


0.38 


16 


1 




7 


287.933 


288.200 


-0.900 


-0.600 


0.4 


0.2 


0.06 


0.48 


0.31 


9 


14 




total 










34.6 


13.2 


5.78 


0.38 


0.44 









^Mi2, Mi3, and Mis represent ^^CO mass, ^^CO mass, and C"*^®© mass, respectively. Units are in lO** Mq. 
^Average value of the T^x of each of the C^^O cores within the region. 



Table 2 

Observed Properties of the C^^O Cores 



Name 


(°) 


Coordinate 
b a 
(°) (J2000) 


S 

(J2000) 


(Ci«0) 
(K) 


Vlsr 
(km s"^) 


AV 
(km s-^) 


/is'' 




Iis/Il3 


1 


285.967 


-0.500 


10 


34 23.6 


-58 


46 


56 


0.6 


-16.6 


2.4 


1.6 


9.5 


0.17 


2 


286.100 


-0.133 


10 


36 41.1 


-58 


31 


50 


0.3 


-20.7 


4.5 


1.4 


10.6 


0.14 


3 


286.100 


0.233 


10 


38 03.5 


-58 


12 


40 


0.5 


-20.6 


2.6 


1.4 


16.5 


0.09 


4 


286.167 


0.033 


10 


37 45.3 


-58 


25 


06 


0.4 


-21.6 


2.6 


1.2 


4.3 


0.28 


5 


286.233 


0.167 


10 


38 41.8 


-58 


20 


04 


1.2 


-20.0 


2.9 


3.6 


17.6 


0.20 


6 


286.400 


-0.233 


10 


38 19.0 


-58 


45 


54 


0.5 


-24.1 


3.1 


1.9 


9.3 


0.20 


7 


286.433 


-0.400 


10 


37 54.6 


-58 


55 


36 


0.6 


-13.0 


2.5 


1.5 


13.4 


0.11 


8 


287.133 


-0.833 


10 


41 03.7 


-59 


38 


42 


0.3 


-19.9 


3.8 


1.2 


12.2 


0.10 


9 


287.233 


-0.533 


10 


42 53.0 


-59 


25 


44 


0.5 


-17.6 


2.6 


1.2 


9.4 


0.13 


10 


287.333 


-0.600 


10 


43 19.7 


-59 


32 


05 


0.3 


-18.7 


4.0 


1.2 


12.7 


0.10 


11 


287.333 


-0.433 


10 


43 56.8 


-59 


23 


16 


0.3 


-19.3 


5.4 


1.5 


17.1 


0.09 


12 


287.700 


-0.733 


10 


45 24.0 


-59 


49 


28 


0.3 


-25.8 


3.9 


1.1 


14.1 


0.08 


13 


287.967 


-1.133 


10 


45 48.2 


-60 


18 


09 


0.7 


-21.3 


2.8 


2.3 


17.2 


0.13 


14 


288.100 


-0.700 


10 


48 21. ,3 


-50 


58 


44 


0.5 


-12.1 


2.2 


1.1 


3.4 


0.32 


15 


288.233 


-1.133 


10 


47 43.0 


-60 


25 


31 


0.4 


-17.8 


3.5 


1.3 


9.7 


0.14 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and 

arcseconds. 

^/i8 and /i3 represent the integrated intensity of the C^^O and ^^CO emission, respectively. Units are in K km s~^. 



Table 3 

Physical Properties of the C^*0 Cores 



Name 


Tex 

(K) 


JV(H2)» 
(10^1 cm-=) 


{N(H2)) 
(10=1 cm- = ) 


r 
(pc) 


'"dec 

(pc) 


M 
(Ms) 


n(H2) 
(cm-=*) 


^ V^omp 

(km s-i) 


(Me) 




IRAS" 


MSX"* 


jjl3cQ+e 


1 


16 


10.6 


7.4 


2.0 


1.8 


2100 


900 


2.7 


3000 


1.4 








2 


15 


9.2 


6.9 


3.1 


2.9 


4600 


540 


3.2 


6500 


1.4 








3 


21 


11.1 


8.4 


2.7 


2.5 


4400 


750 


3.3 


6100 


1.4 








4 


20 


9.1 


7.3 


1.4 


1.0 


1000 


1200 


2.6 


1900 


1.9 








5 


28 


34.8 


25.7 


1.8 


1.6 


6100 


3400 


2.8 


2900 


0.5 


1 


1 


1 


6 


26 


17.3 


12.8 


2.0 


1.8 


3700 


1600 


3.6 


5300 


1.5 


1 


3 


2 


7 


17 


10.7 


8.6 


1.6 


1.3 


1600 


1300 


3.1 


3200 


2.0 








8 


24 


10.1 


7.0 


2.7 


2.5 


3700 


630 


3.6 


7500 


2.0 








9 


21 


9.7 


7.3 


1.6 


1.3 


1400 


1100 


2.5 


2200 


1.6 




1 


3 


10 


25 


11.0 


11.0 


0.8 




520 


3300 


4.0 


2700 


5.2 




2 




11 


22 


11.8 


10.1 


1.4 


1.0 


1400 


1700 


4.5 


6100 


4.3 








12 


20 


8.8 


6.5 


2.3 


2.1 


2600 


710 


3.5 


5900 


2.3 




1 




13 


20 


17.4 


12.8 


2.2 


1.9 


4200 


1400 


4.9 


11000 


2.6 








14 


9 


5.9 


4.8 


1.2 


0.6 


450 


1000 


1.9 


900 


2.0 








15 


19 


9.9 


7.0 


2.0 


1.8 


2000 


850 


2.9 


3500 


1.8 






4 


average 


20 


12.5 


9.6 


1.9 


1.7 


2600 


1400 


3.3 


4600 


2.1 








median 


20 


10.6 


7.4 


2.0 


1.8 


2100 


1100 


3.2 


3500 


1.9 









^The column density for the peak position of each core. 
^Beam-deconvolved radius of each core. 
•^Number of associated IRAS point sources. 
'^Number of associated MSX point sources. 
®The name of the associated H^^CO"*" cores. 



Table 4 

Average Physical Properties of C^®0 Cores in Massive-Star-Forming Regions 



Region 


Distance 


Number of 


Number of 




JV(H2) 


M 




n{^2Y 




Ref. 




(pc) 


C^^O Cores 


Resolved Cores 


(km s"^) 


(10^^ cm-=) 


(Mo) 


(pc) 


(10^ cm-^*) 






Orion B 


400 


19 


18 


1.7 (1.7) 


2.1 (1.6) 


250 (160) 


0.5 (0.5) 


7.6 (5.7) ' 


2.1 (1.4) " 


1 


Orion A 


480 


19 


12 


1.7 (1.6) 


2.4 (1.7) 










2 


Vela C 


700 


13 


10 


2.5 (2.3) 


2.8 (2.8) 










3 


Cepheus OB3 


730 


8 


8 


2.0 (1.9) ^ 


1.5 (1.4) 










4 


S35/37 


1800 


18 


4 


1.7 (1.6) 


1.7 (1.1) 


1100 (640) 


1.3 (1.3) 


11.0 (9.2) 


0.5 (0.6) 


5 


77 Car CMC 


2500 


15 


11 


3.2 (2.9) 


1.3 (1.1) 


2600 (2100) 


2.2 (2.0) 


1.2 (0.9) 


1.6 (1.3) 


6 


Ccntaurus III 


3500 


14 


10 


2.3 (2.3) 


1.3 (1.3) 


5400 (4700) 


2.8 (2.8) 


1.1 (1.0) ° 


0.8 (0.6) 


7 


Ccntaurus II 


4500 


5 


5 


3.1 (3.0) 


1.9 (1.9) 


20000 (24000) 


4.2 (3.8) 


1.1 (0.7) ° 


0.6 (0.5) 


7 


Centaur us I 


5300 


16 


11 


2.3 (2.2) 


1.3 (1.0) 


11000 (8600) 


4.0 (3.6) 


0.7 (0.7) " 


0.5 (0.5) 


7 



Note. — Median values are given in parentheses. 



References. — (1) Aoyama ct al. 2001; (2) Nagahama 1997; (3) Yamaguchi et al. 1999b; (4) Yu et al. 1996; (5) Saito et al. 1999; (6) this work; (7) Saito et 
al. 2001. 

""Number of C^^O cores with rdec/r > 0.8. 
^FWHM line width at the peak position of the core. 
^Values for resolved cores. 

•^Calculated by using the equation (Mvir/M©) = 209 X (r/pc) (AV/km s~^)^. 
^Calculated in this work. 

^Derived from the composite profile, i.e., AV^omp- 



Table 5 

Properties of the Associated IRAS Point Sources 



Coordinate Flux Density 



Name 




b 


a 


S 


Fi2 


^25 






L 


Quality 


c.c.» 






n 


n 


(J2000) 


(J2000) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 










10365-5803 


286.203 


0.170 


10 38 30.6 


-58 19 00 


7.242 


85.64 


1173 


2782 


30000 


1333 


CABA 


5 


10361-5830 


286.375 


-0.255 


10 38 03.8 


-58 46 18 


12.38 


38.38 


625.7 


2156 


21000 


1333 


DCBD 


6 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcminutes, and arcseconds. 

^Correlation coefficient. 

^The name of the associated C""^®© core. 



Table 6 

Properties of the Associated MSX Point Sources 









Coordinate 






Flux Density 








No. 


Name 


I 


b 


a 


S 


Band-A 


Band-C 


Band-D 


Band-E 


Quality 








n 


n 


(J2000) 


(J2000) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 






1 


MSX6C G286.2086+00.1694 


286.2086 


0.1694 


10 38 32.5 


-58 19 12 


1.3532 


2.8816 


7.1817 


40.572 


4444 


5 


2 


MSX6C G286.3579-00.2933 


286.3579 


-0.2933 


10 37 48.4 


-58 47 48 


0.71259 


1.8153 


2.6774 


6.0653 


4344 


6 


3 


MSX6C G286. 3747-00. 2630 


286.3747 


-0.2630 


10 38 02.0 


-58 46 43 


3.5907 


4.7559 


2.4093 


7.5775 


4444 


6 


4 


MSX6C G286. 3773-00. 2563 


286.3773 


-0.2563 


10 38 04.6 


-58 46 26 


1.6278 


2.9183 


3.8554 


12.095 


4444 


6 


5 


MSX6C G287. 2238-00. 5339 


287.2238 


-0.5339 


10 42 48.9 


-59 25 29 


1.2040 


1.6909 


1.5151 


4.3560 


4344 


9 


6 


MSX6C G287. 3399-00. 6008 


287.3399 


-0.6008 


10 43 22.3 


-59 32 19 


0.80028 


2.9386 


5.9723 


12.516 


4444 


10 


7 


MSX6C G287. 3506-00. 5900 


287.3506 


-0.5900 


10 43 29.2 


-59 32 03 


0.55922 


2.0313 


3.9602 


7.9096 


4444 


10 


8 


MSX6C G287. 6393-00. 7219 


287.6393 


-0.7219 


10 45 00.9 


-59 47 10 


3.4863 


4.2647 


3.6021 


17.835 


4444 


12 



Note. — Units of right ascension are hours, minutes, and seconds, and units of decHnation are degrees, arcminutes, and arcseconds. 
^The name of the associated C^^O core. 



Table 7 

Physical Properties of H^^C0+ Cores 



Name 


a 

(J2000) 


5 

(J2000) 


(K) 


Vlsr 
(km s"^) 


Ay 

(km s"^) 


(K) 


r 


JV(H2)* 


r 
(pc) 


^ VI;;omp 

(km s"^) 


M 
(Mo) 


(Mq) 




n(H2)'' 




1 


10 38 31.7 


-58 18 44 


0.38 


-19.4 


2.9 


30 


0.018 


4.5 


0.77 


3.8 


1400 


2300 


3.4 


1.1 


5 


2 


10 37 52.9 


-58 47 20 


0.41 


-22.5 


3.0 


55 


0.011 


7.6 


0.27 


3.0 


400 


500 


7.6 


6.8 


6 


3 


10 42 47.8 


-59 26 24 


0.35 


-17.3 


1.3 


35 


0.017 


2.0 














9 


4 


10 47 46.7 


-60 26 26 


0.81 


-17.8 


1.6 


30 


0.035 


5.0 


0.39 


1.6 


470 


200 


4.4 


2.8 


15 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcminutes, and arcseconds. 

^Units are in 10^^ cm~^. 

''Units are in 10^ cm~^. 

*^The name of the associated C^^O core. 



Table 8 
Properties of the Outflow 

















Kinetic 


Mechanical 




Dynamical 










Size'^ 


Radius 


Mass 


Energy 


Luminosity 


Momentum 


Timeseale 


Component 


(km s 


-h 


(km s"^) 


(pc) 


(pc) 


{Me) 


(10*^ ergs) 


(Le) 


(M0 km s-^) 


(10* yr) 


Blue 


(-35, 


-25) 


16.5 


0.41 


0.41 


1.2 


1.7 


0.39 


14 


2.4 


Red 


(-15, 


-5) 


13.5 


0.41 


0.55 


2.1 


1.9 


0.37 


20 


4.0 



^Velocity range where the high velocity emission is detected at > l-cr noise level (T^-^O.S K). 
^Maximum velocity shift of the """^CO (-7 — 2-1) emission from Vlsr — —18.5 km s"^"'". 
'^The maximum separation of the outflow lobes from the position of the driving source. 



Table 9 

Average Physical Properties of C^*0 Cores in the t] Car CMC 





Group* 


Name*" 


^ VJ;omp 


Ar(H2)'= 


M 




n(H2)'^ 










(km s"^) 




(Mo) 


(pc) 


(x 10^ cm-^*) 




Cores with IRAS point sources 


A 


5, 6 


3.2 (3.2) 


26.1 (26.1) 


4900 (4900) 


1.9 (1.9) 


2.5 (2.5) 


1.0 (1.0) 


Other star-forming cores 


B 


9, (10), 12, 15 


3.2 (3.2) 


9.9 (9.8) 


1600 (1700) 


2.0 (2.0) 


0.9 (0.9) 


1.9 (1.8) 


All star-forming cores 


A-l-B 


5, 6, 9, (10), 12, 15 


3.2 (3.2) 


15.3 (10.5) 


2700 (2300) 


2.0 (2.0) 


1.5 (1.1) 


1.5 (1.6) 


Non-star-forming cores 


C 


1, 2, 3, (4), 7, 8, (11), 13, (14) 


3.3 (3.2) 


10.7 (10.6) 


2600 (2100) 


2.4 (2.4) 


0.9 (0.8) 


1.8 (1.7) 


All C^®0 cores 




all 


3.3 (3.2) 


12.5 (10.6) 


2600 (2100) 


2.2 (2.0) 


1.2 (0.9) 


1.7 (1.6) 



Note. — Median values are given in parentheses. 
^See text (§ 4.1). 

^The Names of the associated C^^O cores. The names of the unresolved cores are given in parentheses. 
•^Units are in 10^^ cm~^. 
"^Values for resolved cores. 



